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Solving dynamical problems in general relativity requires the full machinery of numerical relativity. 
Wilson has proposed a simpler but approximate scheme for systems near equilibrium, like binary 
neutron stars. We test the scheme on isolated, rapidly rotating, relativistic stars. Since these objects 
are in equilibrium, it is crucial that the approximation work well if we are to believe its predictions 
for more complicated systems like binaries. Our results are very encouraging. 
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I. INTRODUCTION 

Some of the most interesting unsolved problems in gen- 
eral relativity require full dynamical solutions of Ein- 
stein's equations in three spatial dimensions. Such so- 
lutions have to be found numerically, and this is only 
barely becoming technically feasible. An important set 
of problems in this category is the binary coalescence of 
black holes and the binary coalescence of neutron stars. 
Such events are expected to be a significant source of 
gravitational waves that will be detectable by new gen- 
erations of detectors such as LIGO. 

In Newtonian physics, binary stars can orbit in an 
equilibrium system. In general relativity, by contrast, 
a binary system loses energy by gravitational wave emis- 
sion. The orbit shrinks, and the two stars ultimately 
coalesce. Though this is clearly not an equilibrium situa- 
tion, the orbital decay occurs on a much longer timescale 
than an orbital period, at least up until the last plunging 
orbit when the stars are very close. Preliminary calcu- 
lations of binary coalescence and gravitational collapse 
suggest that the amount of energy radiated gravitation- 
ally is small. Thus, even when the system becomes highly 
dynamical and far from equilibrium, one might expect 
that it is still the nonradiative part of the gravitational 
field that controls the evolution. 

Wilson [Q,^ has proposed an approximation scheme 
that tracks the evolution of coalescing binary neutron 
stars without solving the full dynamical Einstein field 
equations. The method may also be applicable to binary 
black hole systems The scheme applies to systems 
that are either in or near equilibrium, in which case a 
reduced set of Einstein's equation should adequately de- 
scribe the system. For example, a binary system is near 
equilibrium as long as the emission of gravitational radi- 
ation is small. In strict equilibrium, as in the case of a 



single rotating star, there is a coordinate frame in which 
the first and second time derivatives of the metric are 
zero. In the 3-1-1 formalism, this means in particular 
that the time derivatives of the 3-metric ^ij and the ex- 
trinsic curvature Kij are zero. In quasi-equilibrium, the 
time derivatives are small, and the metric and extrinsic 
curvature will not depart significantly from their initial 
values. Wilson's approximation consists in setting time 
derivatives exactly equal to zero in a selected subset of 
Einstein's equations, and ignoring the remaining dynam- 
ical equations. This approximation results in a smaller, 
more tractable set of field equations. Part of the strategy 
for selecting the subset of Einstein's equations is to guar- 
antee that and Kij are solutions of the initial-value 
(or constraint) equations. Wilson has proposed evolving 
the system through a sequence of initial-value problems 
by solving the full dynamical equations for the matter in 
the instantaneous background metric, and then updating 
the metric quantities at each time step by re-solving the 
selected subset of Einstein's equations. We will outline 
below a simpler method to track the evolution, which 
exploits the near equilibrium of the matter as well. 

As compelling as this idea sounds, it is impossible to 
calibrate the approximation without comparing it with 
solutions to the exact equations. No such exact solutions 
exist for realistic, dynamical 3-dimensional cases. Only 
recently has it become possible to solve Einstein's equa- 
tions numerically for interesting 2-dimensional problems. 
In fact, it is only in the last few years that as simple a 
problem as the equilibrium structure of a rapidly rotat- 
ing relativistic star could be thoroughly investigated. In 
this paper, we use these rotating equilibrium solutions 
to calibrate Wilson's approximation scheme. This is the 
simplest case for which the approximation scheme is dif- 
ferent from the exact equations. Because the system is 
a true equilibrium, it is clearly necessary that the ap- 
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proximation work well in this case. Only then will we 
have confidence that the method is at all useful in more 
complicated situations such as binary systems. 



II. BASIC EQUATIONS 



A general metric may be written in 3 + 1 form as 



'-dt^ + -i,j{dx' + I3'dt){dx^ + 13^ dt). 



The dynamical equation for jij is 



(1) 



(2) 



where Di denotes a covariant derivative with respect to 
7ij . The trace of this equation is 



5t In 7^/2 = ^aK + D,l3' 



(3) 



where 7 = det7y and K = K\. The trace- free part of 
Eq. (D is 

+ D,pj+Djp,^^-,,jDkpK (4) 

We fix the six components of the extrinsic curvature Kij 
by demanding that each data slice be a maximal slice and 
that the left-hand side of Eq. (^) be equal to zero. This 
gives 



and 



if = 0, 



2aK,j = D,I3, + D,I3, - -^..D^pK 



(5) 



(6) 



Note that da ^ unless A/3* = 0. 

To solve the Hamiltonian constraint equation, it is con- 
venient to use a conformal decomposition of the spatial 
metric. To satisfy the demand that the left-hand side of 
Eq. ffl) be zero, we choose the metric to be conformally 
flat so that 7~^/^7ij = /ij, where fij is the flat met- 
ric in whatever coordinate system is used. Therefore, we 
decompose the spatial metric as 



7»j 



(7) 



The conformal factor $ is determined then by the Hamil- 
tonian constraint 
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where the source term is 



(8) 



(9) 



Here n" is the normal vector to a i = constant slice. 
Tab is the stress-energy tensor, and is the flat-space 



Laplacian. Note that although indices range over 

1, . . . , 3, indices a, 6, . . . range over 0, . . . , 3. 

The shift vector is determined by substituting Eq. (||) 
into the momentum constraint 



where 



be 



(10) 



(11) 



We use the results that for a conformally flat metric we 
may write 



and for if = 0, 



, (12) 



(13) 



where Vj denotes the covariant derivative in flat space. 
Thus Eq. becomes 



(14) 



V2/3* + xV'(Vj/3^') = 16™$^ 5* 



This equation can be simplified to two equations, one in- 
volving a vector Laplacian and the other a scalar Lapla- 
cian, by setting 



(3' = G' - \V'B. 
4 



(15) 



Then the two equations that must be solved are 
V^G* = 16Tra<i>^S' (16) 



and 



(17) 



Though we are not imposing the full set of dynami- 
cal equations for the evolution of Kij, we do have the 
freedom to preserve the maximal slicing condition by 
requiring dtK = 0. The resulting equation can also be 
written with a simple Laplacian by using Eq. (H). The 
result is the lapse equation 



V2(a$) = (a$) 
where 



-<^^K,jK'^ + 27r$'*(/9 + 25) 



S = Y'T,, 



(18) 



(19) 
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The above field equations, in combination with the 
matter equations to be discussed below, form a cou- 
pled nonlinear set that must be solved by iteration. The 
boundary conditions for the field quantities follow from 
asymptotic flatness; the speciflc form depends on the ap- 
plication. We are especially interested in uniformly ro- 
tating configurations such as binary neutron stars in syn- 
chronous orbit. For such systems we work in a corotating 
coordinate system so that there is no time variation of 
the fields (in the near equilibrium approximation of the 
method) . Following Wilson j2| , we can implement this 
by replacing Eq. (15) with 
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^ G - -\I'B + (J^ X r) 



(20) 



which leaves Eqs. (17) and ( [16D unchanged. Here $7 is 
the constant angular velocity of the system. 

For the matter, we will consider a perfect fluid for 
which 



Tab = (PO + P^ + P) UaUb + Pgab 



(21) 



Here po is the rest-mass density, pi is the internal energy 
density, P is the pressure, and J/" is the fluid 4- velocity. 
For this source, the density p in Eq. (||) is 

P = (po + P^ + P) {aU')^ - P, (22) 

the momentum source in Eq. ( pi] ) is 

^{po + p^ + P){aU')f^U„ (23) 

and the source term 5* in Eq. ( p^ is 

S^ipo + A + P) [(«C/*)' - 1] + 3P. (24) 

We treat fluids that are in uniform rotation, for which 
the 4-velocity J/" is given by 



dt d(f) 



The normalization condition U ■ U = — 1 gives 



1 



1/2 



(26) 



Now consider the equations for the matter in the near 
equilibrium approximation. The key approximation is 
that in the corotating frame there is a Killing vector that 
is timelike everywhere. In the nonrotating coordinates, 
this vector can be written as 



(27) 



Because the 4-velocity (^5|) is proportional to a Killing 
vector, the matter equations may be integrated to give 
the hydrostatic equilibrium result 



— = constant, 
h 



where 



hih 



dP 



(28) 



(29) 



Po +Pi + P 
For a polytropic equation of state 

P = Kp^o, (30) 
where K and F are constants, we have 



P 



Pi = 



F- 1 



h = 



Po + Pi+ P 
Po 



(31) 



In this approximation, we have reduced all of the hydro- 
dynamics to a single algebraic equation, Eq. (p8|). 



III. AXISYMMETRIC ROTATING STAR: 
EQUATIONS 

To calibrate the method, we apply it to a true equi- 
librium system in axisymmetry and compare with the 
complete numerical solution found with no approxima- 
tions. For this purpose, we use models of rotating neu- 
tron stars supported by a polytropic equation of state. 
Fully relativistic models have been constructed by several 
authors (see Refs. |^-^ and references therein). Solving 
Einstein's equations for these stars is nontrivial numeri- 
cally. It is only the recent availability of such solutions 
that makes this calibration feasible. 

In spherical polar coordinates and axisymmetry, we 
find that Eqs. (|l6|) and ( [T7| ) are satisfied by setting the 
quantity B of Eq. (|l^) to zero and with the only nonzero 
component of the shift vector f3'^ = (3. Note that this im- 
plies, not only that the left-hand side of Eq. (^) is zero, 
but also that (9t7 = 0. This means that we are find- 
ing a stationary solution of the approximate equations. 
Given this solution for the shift vector, the term K^K"^^ 



(25) appearing in Eqs. (||) and ([T^) is given by 



sm 



2a2 



(32) 



where commas denote partial derivatives. Only the 0- 
component of the vector Eq. (|l^) is nontrivial, and be- 
comes the scalar equation 



2_a 

r dr 



2 cot 6* d 



r2 86 
1 da 



(3 



a dr ^ dr J dr 



6 8(3 167ra 



ade ^ 86 ) de 



r-2 sin^ 9 



S^. (33) 



The 4-velocity components appearing in the matter 
sources are given by 



t/* = [^2 _ sin^ 6{(3 + nf] 



21-1/2 



(34) 
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The above equations turn out to be simplified versions 
of the exact equations for stationary, axisymmetric con- 
figurations given in Ref. § (henceforth CST §). The 
exact metric has four nonzero metric coefficients, denoted 
by 7, p, a and lu in CST, though the approximate met- 
ric here has only three: a, /3, and Thus even though 
there is no dynamics in the field, and even though the 
equation of hydrostatic equilibrium for the matter is rig- 
orously obeyed, the Wilson scheme is still an approxima- 
tion for this problem. The correspondence between the 
approximate and exact metric coefficients is given by 



(35) 



The fiuid velocity v in the ZAMO frame used in CST is 
given by 



1 



1 



(36) 



In spherical symmetry, the approximate scheme reduces 
to the exact scheme, with two nonzero metric coefficients. 
We will now quantify the degree of error in the nonspher- 
ical axisymmetric case. 

We can take over the numerical scheme of CST to 
solve the approximate equilibrium equations. In fact, 
the structure of the equations is very close in that they 
involve the same differential operators on the left-hand 
sides. In particular, Eqs. (||) and ( [l8| ) involve V^, as 
in Eq. (3) of CST, and the operator in Eq. (|3^) is the 
same as that in Eq. (5) of CST. Therefore the solution is 
computed as in Eqs. (27) and (29) of CST. The nondi- 
mensionalized source terms analogous to Eq. (30) of CST 
are 



<S'$(s,/i) 



1 <I>^ 

16 (a$)^ 

X |[(1 ^ *) so), 

2 



(37) 



5^2 



27r$''f: 



S'q$(s, /i) = 



4^2 



I- S 



7 $^ 



Po + Pi + P) 



p 



1 - s 



16 (a$)2 
x{[(l-,s).ci,,]^ + (1-^*2)^;^^} 

1 



(38) 



+ 47r$ r 



1 - s 



1 - s 



{po +P^+P)Y 



Po + P^ + P) 



P 



3P 



and the source term analogous to Eq. (32) of CST is 



1 , X 7 



(39) 



Po + p» -f P 

l-s) l-v^ 



where s is an auxiliary radial coordinate defined in CST. 
The entire iterative scheme used to solve the approximate 
equations is identical to the one in CST. 

To calibrate the approximation, we first compute an 
exact sequence of constant rest mass polytropes of in- 
creasing angular momentum. Each member of the se- 
quence is specified by two parameters: the ratio of polar 
to equatorial radius, and the central rest-mass density. 
We next compute the approximate sequence using the 
same values for these two parameters for each model. 
We then compare the metric coefficients of correspond- 
ing models using the relationships in (|35|). We also com- 
pare global quantities such as the total mass and angular 
momentum. As a further diagnostic, we calculate two 
relativistic virial quantities |l0|,|ll[ whose values should 
be identically one for an exact equilibrium solution. In 
the notation of CST, these quantities are 



Aorf = 327r 



P + (e + P) 



e^°'rdrd0 



(40) 




dj dp\^ ^ 1 f ^ dp^ ^ 
dr dr 



3e-2''sin2( 



2 



dr J \ dtheta J 



3P + (e + P) 



1 -t>2_ 



>rdrd6. 



(41) 



d{l + p)9(7 + p) — dad"f + dadp 



da 



57 



r tan 9 89 partialr 
1 ^ dp ^ 1 dp 



r tan 9 09 dr r tan 9 89 



-^Pr^ sin^ 9dujduj 



= (7-p)/2^2 



where 



dadp 



da dp ^ 1 da dp 
dr dr r'^ dr d9 



r sin 9 drd9, 



(42) 



and e ^ Po + Pi is the total mass-energy density. Here 
X^d involves an integration with a 3-dimensional volume 
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element sin 9 dr d9 and is the relativistic generalization 
of the classical virial theorem 



2^kin + 3(r - l)^int + C/grav = 0. 



(43) 



The quantity \2d involves an integration with a 2- 
dimensional volume element r dr dO. The discrepancy 
from unity is a measure of numerical error for our so- 
lutions of the exact equations. It is a measure of the 
larger inherent error for our solutions of the approximate 
equations. 



IV. AXISYMMETRIC ROTATING STAR: 
NUMERICAL RESULTS 

To calibrate the approximate scheme against the ex- 
act solution, we choose the most stringent case, in which 
the configuration is very relativistic and rapidly rotating. 
When it is rotating rapidly, there are large deviations 
from spherical symmetry, so that the approximation is 
no longer exact. For polytropes, the largest rotation is 
attained for nearly incompressible matter, i.e. for large 
r = 1 + 1/ri or small polytropic index n. We choose 
n = 0.5. 

In constructing an exact sequence of rotating equilib- 
ria as a benchmark, we start with a nonrotating star 
having a central value of energy density e = 1 (note 
that all "barred" quantities are nondimensional as de- 
fined in CST). This configuration is relativistic, with 
M/R = 0.298 and rest mass Mq = 0.148, just below 
the maximum rest mass of a nonrotating star for this 
equation of state (Mq = 0.151). Holding the rest mass 
constant, we construct a sequence of increasing uniform 
rotation, up to the point of mass shedding. As described 
above, we then construct the corresponding models with 
the same central value of e and ratio of polar to equato- 
rial radius using the approximate scheme. A comparison 
of some of the global quantities for the sequence is given 
in Table |. The high values of polar redshift Zp and T /W 
confirm that the sequence is both highly relativistic and 
rapidly rotating. As expected, the deviations are largest 
near the mass shed limit, but even there they are never 
worse than about 1%. 

We can understand why the overall discrepancy is 
small by looking at Fig. |l|. Here we plot a measure of 
the deviation in the exact solution from conformal flat- 
ness, which is assumed in the approximate method. In 
the figure we plot the angular profile at selected radii of 
the quantity 



A 



_ acsT - (7 - /3)/2 



acsT 



(44) 



computed for the exact rotating model with T/W = 
0.159. Note that this quantity is identically zero on the 
axis because of local fiatness there. The maximum devi- 
ation occurs on the equator (f = 0.48), but is only about 
5%. Outside the star, A — > as r — > 00. 
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FIG. 1. Angular profile of the deviation of the exact solu- 
tion from conformal flatness at selected radii. The deviation 
A is defined in Eq. (^). The star is a rapidly rotating, highly 
relativistic polytrope with n — 0.5 and rest mass just below 
the maximum rest mass of a nonrotating star for this equa- 
tion of state. The radii f are in the nondimensional units of 
CST, and = cos6l. 

In Fig. H we plot along an equatorial radius the frac- 
tional error in the conformal factor. 



(5$ 



(45) 



where $oxact = exp[(7 — p)/4]. Similarly, in Fig. ||we plot 
the fractional error Suj. Figure Q shows the mass-energy 
e along an equatorial radius for the two schemes. The 
two profiles are almost coincident. 

A further comparison is provided by Fig. |^, which 
shows the virial quantities X2d and X^d along the se- 
quence, computed for each of the two schemes. In the 
case of the exact method, the deviation from unity is 
a measure of numerical error, which is less than 0.1%. 
The deviation for the approximate method measures the 
inherent error, which is about a factor of 10 bigger. 

To push the approximate scheme to the limit, we 
now consider a second equilibrium sequence, a "supra- 
massive" sequence. This sequence has no nonrotating 
member, since its rest mass exceeds the maximum rest 
mass of a nonrotating star for this equation of state 
(Mq = 0.151). Thus the sequence exists only by virtue of 
rotation. We construct the supramassive sequence with 
Mq = 0.176. We expect the discrepancy between the 
approximate and exact methods to be somewhat larger 
for this sequence since it is everywhere far from spherical 
symmetry. This expectation is borne out in Table || and 
Fig. ||. Nevertheless, the discrepancy is not very large. 
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FIG. 2. Fractional error in the conformal factor "I> along an 
equatorial radius for the star in Fig. |ll 
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FIG. 4. Total mass-energy density e along an equatorial 
radius for the star in Fig. gj. The solid line shows the exact 
solution, the dotted line, the approximate solution. 
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FIG. 3. Fractional error in the metric coefficient uj along 
an equatorial radius for the star in Fig. ll. 
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FIG. 5. Virial quantities along the sequence in Table ^ 
Results for the exact equations are shown by the solid line for 
\2d and the dotted line for X^d- Results for the approximation 
are shown by the short dash line for A2d and the long dash 
line for Asd. 
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FIG. 6. Virial quantities along the supramassive sequence 
in Table ^ Results for the exact equations are shown by the 
solid line for \2d and the dotted line for Asd. Results for the 
approximation are shown by the short dash line for \2d and 
the long dash line for X^d- 



V. CONCLUSION 

We have tested Wilson's approximation scheme on 
rapidly rotating relativistic stars. Since these are equilib- 
rium objects, it is necessary that the scheme give reason- 
ably accurate results if we are to believe its predictions 
for more complicated systems such as binaries. In fact, 
we have found that the method works remarkably well, 
even for highly relativistic objects far from spherical sym- 
metry. The largest errors in any quantities we examined 
were around 5%, and in general the errors were much 
smaller. Global measures such as virial quantities were 
in error by far less than 1%. This agreement is very en- 
couraging. 
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TABLE I. Quantities characterizing a "normal" evolutionary sequence of n=0.5 polytropic neu- 
tron star models. For each value of the central energy density e we first display the results from 
solving the exact equations. Below this, indicated with a dash in the energy density column, are 
the results obtained by solving the approximate equations. 
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0, 


,1132 


0, 


,5252 


0, 


.712 


0, 


.6616 





.4076 


0, 


,5527 


0.5417 


0, 


.7224 


0, 


.01680 


0. 


,1287 


0, 


,1484 


0, 


,1416 


0, 


.5743 


0, 


.778 


0, 


.6440 


0, 


.4552 


0. 


,5422 







.7256 





.01665 


0. 


,1288 


0, 


,1485 


0, 


,1413 


0, 


.5710 





.782 





.6445 





.4550 


0, 


,5426 


0.5148 





.7226 





.01795 


0. 


,1293 


0, 


,1484 


0, 


,1559 


0, 


,6124 





.817 





.6349 


0, 


.4875 


0, 


,5348 




0, 


.7260 


0, 


.01777 


0. 


1293 


0, 


,1485 


0, 


,1555 


0, 


.6091 


0, 


.820 


0, 


.6354 


0, 


.4876 


0, 


,5351 


0.5115 


0, 


.7221 


0, 


.01812 


0. 


,1294 


0. 


,1484 


0. 


,1577 


0, 


.6217 


0, 


.824 


0, 


.6337 


0, 


.4956 


0. 


,5337 




0, 


.7254 


0, 


.01792 


0. 


,1294 


0. 


,1485 


0. 


,1572 


0, 


.6184 


0, 


.827 


0, 


.6341 


0, 


.4958 


0. 


,5339 


0.5098 





.7218 





.01821 


0. 


,1294 


0, 


,1484 


0, 


,1587 


0, 


.6303 





.830 





.6331 





.5034 


0, 


,5332 







.7249 


0, 


.01799 


0. 


,1294 


0, 


,1484 


0, 


,1581 


0, 


,6272 





.833 


0, 


.6334 





.5037 


0, 


,5331 


0.5094 


0, 


.7216 





.01822 


0. 


,1294 


0, 


,1484 


0, 


,1588 


0, 


.6331 


0, 


.832 


0, 


.6329 


0, 


.5060 


0, 


,5330 







.7248 





.01800 


0. 


,1294 


0, 


,1484 


0, 


,1581 


0, 


.6300 


0, 


.835 





.6332 





.5063 


0, 


,5328 


0.5094 





.7216 





.01822 


0. 


,1294 


0, 


,1484 


0, 


,1589 


0, 


,6332 





.832 





.6329 





.5061 


0, 


,5330 







.7248 





.01800 


0. 


,1294 


0, 


,1484 


0, 


,1581 


0, 


,6302 





.835 





.6332 





.5065 


0, 


,5328 


0.5094 


0, 


.7216 


0, 


.01822 


0. 


,1294 


0, 


,1484 


0, 


,1589 


0, 


,6333 


0, 


.832 


0, 


.6329 


0, 


.5062 


0, 


,5330 




0, 


.7248 


0, 


.01800 


0. 


,1294 


0, 


,1484 


0, 


,1581 


0, 


.6302 


0, 


.835 


0, 


.6332 


0, 


.5066 


0, 


,5328 



'^Central energy density. 

''Angular velocity measured at infinity. 

'^Moment of inertia. 

Total mass- energy. 
■^Rest mass. 

'Rotational kinetic energy over gravitational binding energy. 
* Circumferential radius. 

''Eccentricity. 

'Measure of frame dragging. 
■'Matter velocity at equator. 
''Polar redshift. 
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TABLE II. Quantities characterizing a "supramassive" evolutionary sequence of n=0.5 polytropic 
neutron star models. Entries are as described for Table ffl. 



— a 
Ec 




it 




rc 

1 




M 










h 

e 










V k 

Zjp 


1.0957 


0, 


.9464 


0, 


.01552 


0. 


1471 


0. 


,1758 


0. 


,1248 


0, 


.4819 


0, 


.703 


0, 


.8344 


0, 


.5196 


0, 


.9282 




0, 


.9549 


0, 


.01525 


0. 


1472 


0. 


,1760 


0. 


,1246 


0, 


.4760 


0, 


.714 


0, 


.8355 


0, 


.5192 


0, 


.9294 


1.0602 


0, 


.9362 


0, 


.01571 


0. 


1472 


0. 


,1758 


0. 


,1254 


0, 


.4857 


0, 


.705 


0, 


.8284 


0, 


.5174 


0, 


.9132 




0, 


.9444 


0, 


.01544 


0. 


1472 


0. 


,1760 


0. 


,1251 


0, 


.4799 


0, 


.717 


0, 


.8295 


0, 


.5170 


0, 


.9144 


1.0258 


0, 


.9269 


0, 


.01592 


0. 


1472 


0. 


,1758 


0. 


,1264 


0, 


.4901 


0, 


.709 


0, 


.8223 


0, 


.5163 


0, 


.8989 




0, 


.9349 


0, 


.01566 


0. 


1473 


0. 


,1760 


0. 


,1261 


0, 


.4843 





.720 


0, 


.8234 


0, 


.5158 


0, 


.9001 


0.9925 


0, 


.9182 


0, 


.01615 


0. 


1473 


0. 


,1758 


0. 


,1277 


0, 


.4948 


0, 


.714 


0, 


.8162 


0, 


.5158 


0, 


.8851 




0, 


.9260 


0, 


.01588 


0. 


1473 


0. 


,1760 


0. 


,1274 


0, 


.4892 


0, 


.724 


0, 


.8172 


0, 


.5153 


0, 


.8861 


0.9603 


0, 


.9105 


0, 


.01640 


0. 


1474 


0. 


,1758 


0. 


,1294 


0, 


.5003 


0, 


.720 


0, 


.8100 


0, 


.5166 


0, 


.8719 




0, 


.9182 


0, 


.01614 


0. 


1474 


0. 


,1760 


0. 


,1291 


0, 


.4946 


0, 


.730 


0, 


.8110 


0, 


.5162 


0, 


.8730 


0.9292 


0, 


.9031 


0, 


.01668 


0. 


1475 


0. 


,1758 


0. 


,1314 


0, 


.5062 


0, 


.726 


0, 


.8037 


0, 


.5179 


0, 


.8591 




0, 


.9107 


0, 


.01641 


0. 


1476 


0. 


,1760 


0. 


,1311 


0, 


.5005 





.736 





.8048 


0, 


.5176 


0, 


.8602 


0.8991 


0, 


.8964 


0, 


.01698 


0. 


1476 


0. 


,1758 


0. 


,1337 


0, 


.5127 


0, 


.734 


0, 


.7974 


0, 


.5203 


0, 


.8467 




0, 


.9038 


0, 


.01671 


0. 


1477 


0. 


,1760 


0. 


,1334 


0, 


.5071 


0, 


.743 


0, 


.7985 


0, 


.5200 


0, 


.8478 


0.8699 


0, 


.8901 


0, 


.01731 


0. 


1478 


0. 


,1758 


0. 


,1365 


0, 


.5202 


0, 


.742 


0, 


.7911 


0, 


.5239 


0, 


.8348 




0, 


.8975 


0, 


.01703 


0. 


1479 


0. 


,1760 


0. 


,1362 


0, 


.5145 


0, 


.751 


0, 


.7922 


0, 


.5237 


0, 


.8361 


0.8417 


0, 


.8839 


0, 


.01766 


0. 


1480 


0. 


,1758 


0. 


,1395 


0, 


.5284 


0, 


.752 


0, 


.7847 


0, 


.5283 


0, 


.8231 




0, 


.8912 


0, 


.01738 


0. 


1481 


0. 


,1760 


0. 


,1392 


0, 


.5228 





.760 





.7858 


0, 


.5282 


0, 


.8244 


0.8144 


0, 


.8779 


0, 


.01805 


0. 


1482 


0. 


,1758 


0. 


,1429 


0, 


.5380 


0, 


.763 


0, 


.7783 


0, 


.5343 


0, 


.8117 




0, 


.8852 


0, 


.01776 


0. 


,1483 


0. 


,1760 


0. 


,1426 


0, 


.5324 


0, 


.771 


0, 


.7794 


0, 


.5343 


0, 


.8129 


0.7880 


0, 


.8719 


0, 


.01847 


0. 


,1484 


0. 


,1758 


0. 


,1466 


0, 


.5493 


0, 


.775 


0, 


.7719 


0, 


.5423 


0, 


.8004 




0, 


.8790 


0, 


.01817 


0. 


,1485 


0. 


,1759 


0. 


,1462 


0, 


.5438 


0, 


.782 


0, 


.7729 


0, 


.5424 


0, 


.8015 


0.7625 


0, 


.8654 


0, 


.01894 


0. 


,1487 


0. 


,1758 


0. 


,1506 


0, 


.5645 


0, 


.790 


0, 


.7654 


0, 


.5546 


0, 


.7891 




0, 


.8723 


0, 


.01862 


0. 


,1487 


0. 


,1759 


0. 


,1500 


0, 


.5591 


0, 


.797 


0, 


.7663 


0, 


.5549 


0, 


.7898 


0.7593 


0, 


.8648 


0, 


.01901 


0. 


,1487 


0. 


,1758 


0. 


,1512 


0, 


.5677 


0, 


.793 


0, 


.7646 


0, 


.5579 


0, 


.7878 




0, 


.8716 


0, 


.01868 


0. 


,1487 


0. 


,1759 


0. 


,1506 


0, 


.5623 


0, 


.800 


0, 


.7655 


0, 


.5583 


0, 


.7884 


0.7562 


0, 


.8637 


0, 


.01907 


0. 


,1487 


0. 


,1758 


0. 


,1517 


0, 


.5701 


0, 


.795 


0, 


.7637 


0, 


.5598 


0, 


.7863 




0, 


.8705 


0, 


.01873 


0. 


,1487 


0. 


,1759 


0. 


,1510 


0, 


.5647 


0, 


.802 


0, 


.7646 


0, 


.5602 


0, 


.7867 


0.7531 


0, 


.8629 


0, 


.01914 


0. 


,1488 


0. 


,1758 


0. 


,1522 


0, 


.5740 


0, 


.799 


0, 


.7629 


0, 


.5639 


0, 


.7849 




0, 


.8695 


0, 


.01878 


0. 


,1487 


0. 


,1758 


0. 


,1514 


0, 


.5687 


0, 


.805 


0, 


.7637 


0, 


.5644 


0, 


.7851 


0.7500 


0, 


.8620 


0, 


.01921 


0. 


,1488 


0. 


,1758 


0. 


,1528 


0, 


.5818 


0, 


.806 


0, 


.7621 


0, 


.5729 


0, 


.7835 




0, 


.8682 


0, 


.01881 


0. 


,1487 


0. 


,1757 


0. 


,1517 


0, 


.5767 


0, 


.812 


0, 


.7627 


0, 


.5736 


0, 


.7828 



^Central energy density. 

'^Angular velocity measured at infinity. 

'^Moment of inertia. 

Total mass- energy. 
"^Rest mass. 

'Rotational kinetic energy over gravitational binding energy. 

^Circumferential radius. 

''Eccentricity. 

'Measure of frame dragging. 
■■Matter velocity at equator. 
"^Polar redshift. 
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